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Abstract We consider the eigenvalue equation for the largest eigenvalue of certain kinds of non-compact 
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1. Introduction 

This article is concerned with eigenvalue equations on L 1 (J) of the form 

r{x)p[x) + J [u(x,y)p(y) - u(y,x)p(x)] dy = Xp(x) for all x £ I. (1-1) 

Here, p is a probability density on the set /, which is either taken to be a compact 
interval [a, b] or the real line R, i.e., p € L 1 (7) with p > and Jjp(x)dx = 1. Sufficient 
conditions for the existence and uniqueness of solutions of (1.1) were given by Burger, 
see [4, Ch. IV. 3], in which case A is the largest eigenvalue. 

If one is interested in a discrete approximation of (1.1), one faces the problem that the 
operator acting on p is the sum of a multiplication operator and a kernel operator; and 
the former is never compact (apart from trivial cases). Therefore, a direct application of 
most standard methods of approximation theory fails because, for these, compactness is a 
prerequisite. In this article, it will be shown that, under some moderate extra conditions, 
these methods can nevertheless be applied. 

One motivation to study equations of the form (1.1) is their occurrence in population 
genetics, which is concerned with the (micro)evolution of the genetic composition of pop- 
ulations. For many situations, individuals are adequately described by a continuous scalar 
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variable, representing, for example, a quantitative character under selection. This leads 
to the definition of so-called continuum-of-alleles (COA) models, in which individuals are 
identified with this variable, referred to as their type. Usually, selection is then modeled 
by type-dependent fitness values, whereas mutation is described, for every source type, 
by a probability distribution for the mutant types. For a recent review of and relevant 
literature on COA models, see [4]. 

In population genetics, evolution may quite generally be assumed to proceed in contin- 
uous time, with overlapping generations, or in discrete generations. For the COA model, 
in both cases, equilibrium is described by an equation of the form (1.1), compare [12]. 
Here, / is the set of possible types. Assuming the population to be effectively infinite, we 
represent it by the probability density p. 

The notation chosen here best fits the case of continuous time, where r(x) describes 
the effective reproduction rate of type x (i.e., the difference of its birth and death rate), 
the so-called Malthusian fitness, and u(x,y) is the mutation rate u\{y) of type y times 
the density m(x,y) of mutant types x, conditioned on a mutation to occur for y. With 
discrete generations, r(x) has the interpretation of the expected number of offspring of 
an individual of type x, i.e., its Wrightian fitness, and mutation is assumed to occur 
during reproduction with some probability /u(y) for type y. The distribution of mutant 
types is again given by m(x, y), hence u(x, y) — m(x, y) n{y) r(y). In both cases, A equals 
the equilibrium mean fitness J f r(x)p(x) dx. 

There are several reasons why it is desirable to approximate a COA model by a model 
with discrete types. One reason is the need for numerical investigations of COA models, 
since most of them arc not tractable analytically. These inevitably require a discrete 
formulation of the model. Another reason is that recently a simple characterization of the 
equilibrium of discrete mutation-selection models has been found [8] (see also [7, 6, 2]); 
this takes the form of a scalar maximum principle in a limit of infinitely many types that 
densely fill a compact interval. Gaining a better understanding of the relation between 
models with discrete and continuous types is therefore promising to enable a transfer of 
some of these results. 

This article starts with a summary of Burger's results on (1.1) in Section 2, since these 
form the basis for our treatment. We will then consider two methods to approximate 
compact kernel operators and extend them to our case. One, the Nystrom method, is 
applicable to continuous functions r and u on compact intervals I and involves sampling* 
of these functions. This is presented in Section 3. The other one, the Galerkin method, 
is based on projections to finite-dimensional subspaces and works — in principle — for a 
broad class of compact operators. In our case, however, one has to make relatively strong 
assumptions, e.g., that the functions r and u are, in some sense, uniformly continuous. 
Then, it turns out that the local averaging in the projection process can be replaced by 
sampling again (if an additional condition is satisfied). This is discussed in Section 4. 
A comparison of both methods in Section 5 and an outlook in Section 6 complete this 
article. 

* The term sampling is used in the meaning also used in signal processing: Instead of a continuous 
function one considers its values at a (properly chosen) finite set of points. 
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2. General properties 

Let us first put the equilibrium condition (1.1) in operator notation. Since we are inter- 
ested in probability densities, we will consider L 1 (7), or a subspace thereof, as the under- 
lying function space. We define the total mutation rate of type x as 

ui{x) = j u{y,x)dy (2.1) 

and, for notational brevity, 

w = Mi — r . 

Then, (1.1) is equivalent to the eigenvalue equation 

(A + \)p = 0, (2.2) 

where, for elements / of the function space and all x E I, 

(Tf)(x)=w(x)f(x), (2.3) 

{Uf)(x) = Ju(x,y)f(y)dy, (2.4) 

A = T-U. (2.5) 

As mentioned above, being a (non-zero) multiplication operator, T cannot be compact 
(compare [16, Thm. 2.1]). Strong results like analogs to the Perron-Frobenius theorem, 
however, are only available for compact, or at least power compact*, operators, see Schae- 
fer [20, Ch. V]. Therefore one considers the following family of kernel operators: 



(K a f)(x) = J k a (x,y)f{y)dy, 



where 

k a (x,y) - 



w(y) + a ' 

These are, under conditions that will be given shortly, power compact or even compact. 
Their connection to the operator A from (2.5) is stated in the following 

Lemma 2.1. [3, Prop. 2.1(i)] Let T, U be operators in a Banach space X, with U 
being bounded, T densely defined, i.e., D(T) = X, and T + a invertible. Then f is an 
eigenvector of A = T — U with eigenvalue —a, i.e., D{A) = D(T) and 

(A + a)f = 0, 

if and only if g = (T + a)f is an eigenvector of K a = U(T + a)^ 1 with eigenvalue 1, 

(K a - l)g = . 

* An operator is said to be power compact if one of its powers is compact. 
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So, explicitly in our case, the eigenvalue equation (2.2) is equivalent to 



(K x -l)q = 



(2.6) 



with q = (T + X)p. This equation can now be used to find sufficient conditions for the 
existence and uniqueness of a solution of (2.2). 

An important class of bounded kernel operators from L 9 (/) into L p (/) (1 < p, q < oo) 
are the Hille-Tamarkin operators, see [11, Sec. 11.3]. Their kernels need to satisfy 



where (Kf)(x) = f I k(x,y)f(y)dy, k(x,.) denotes the function y \— ► k(x,y), and q' is 
the conjugate exponent to q satisfying ^ + ^- = 1, 1 < q' < oo. The Hillc Tamarkin 
norm \.\ pq turns the set H pq (I) of all Hillc Tamarkin operators into a Banach space [11, 
Thm. 11.5]. Here, we are interested in p = q = 1, in which case (2.7) yields 



and K 2 is compact for every K e Hu(I) [11, Thm. 11.9]. 

Let us now turn to kernel operators that are power compact, positive, and irreducible. 
An operator is called positive if it maps the set of non-negative functions into itself, 
for which, in the case of kernel operators, non-negativity of the kernel is necessary and 
sufficient [11, p. 122]. A kernel operator is irreducible if its kernel satisfies [20, Exm. 4 



/ / k(x, y) Axdy > for all measurable J C I with J|, \I\J\ > 0. 
Ji\j J J 

Here, |J| denotes the Lebesgue measure of a measurable set J. Then, the theorem of 
Jcntzsch [20, Thm. V.6.6], which parallels the Perron Frobcnius theorem for matrices, 
states that the spectral radius is an algebraically simple eigenvalue with an (up to normal- 
ization) unique positive eigenfunction (i.e., strictly positive a.e.*) and the only eigenvalue 
with a positive eigenfunction. 

In our case, the following requirements are sufficient for the K a to be Hille-Tamarkin 
operators [3, Sec. 3]. 

(Ul) u is non-negative and measurable. 

(U2) u\(x) from (2.1) exists for a.e. iel and u\ G L°°(I), i.e., u\ is essentially bounded. 
(By Holder's inequality, this implies that U is bounded, cf. [3, Prop. 3.1(h)].) 

(Tl) w = ui — r is measurable and satisfies ess inf xe i w(x) = 0. (The latter can be 
achieved, without loss of generality, by adding a suitable constant to r.) 



\ K \ P q : = ll fc illp < 00 witn h(x) = \\k(x,.)\\ 



(2.7) 



1-ftTln = / ess sup \k(x, y)\ dx < oo 
Ji yei 



in Sec. V.6] 




* The abbreviation 'a.e.' stands for 'almost every' or 'almost everywhere' and means that the set at 
which the condition it refers to is not fulfilled has zero (Lebesgue) measure. 
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(T2) (w + l)- 1 £ L°°(7) is then already a consequence of (Tl). 

(U4) Jj esssup ygJ u(x, y)/(w{y) + a) dx < oo for one (and then for all) a > 0. 
For a > 0, K a is irreducible if U is [3, proof of Thm. 2.2(c)], i.e., 

/ / u(x,y)dxdy > for all measurable J C I with |J|, \I\J\ > 0. (2.8) 
Ji\j J j 

To keep the equilibrium distribution from having atoms, we assume that there is a set 
J C I with positive measure for which essinf^gj w(x) = such that 

cssinf u(x,y) / (w(x))^ 1 dx > 1 (2.9) 
x,ye.J J j 

or the integral diverges [4, cond. 3" in Sec. IV. 3]. 
Putting everything together, we have the following 

Theorem 2.2 (Burger). Under the above conditions, (1.1) has a unique positive 
solution p £ L 1 (i) with \\p\\i = 1, for which A > is the largest spectral value of —A 
from (2.5). 

Proof. See the above, [3, Thm. 3.5], and [4, Sec. IV.3]. □ 

Note that, due to (Tl), p is positive if and only if q = (w + ot)p is, for a > 0. 
Another result that will be needed in the sequel is 

Lemma 2.3. [3, Lemmas 1-3 and Thm. 2.2(h)] Under the above conditions, the 
spectral radius p(K a ) is, as a function of a, strictly decreasing and satisfies p(K\) = 1 
as well as lim^oo p(K a ) — 0. Thus, p(K a ) < 1 implies a > A and p{K a ) > 1 implies 
a < A. 

Throughout the rest of this article, all the above criteria are assumed to be satisfied, 
namely (Ul), (U2), (U4), (Tl), (T2), (2.8), and (2.9). 

3. Discretization — compact interval 

Let the interval I be compact and C(J) denote the Banach space of bounded, contin- 
uous functions equipped with the supremum norm WfW^ = sup xeI We consider 
operators K of the form 

(Kf)(x) = Jk{x,y)f{y)dy for all x £ I (3.1) 

with a continuous kernel k : I x / — > R. First note these two basic results: 
Proposition 3.1. Any K of the form (3.1) maps L 1 ^) into C(J) C L 1 ^). 
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Proof. Wc follow the proof of [5, Thm. 2.1], where this is shown for L 2 (7), which, 
since I is compact, is a subspace of L 1 (J). Let / e L 1 (7) and x, £ E I be given. Then, 

\{Kf){x) - (Kf)(0\ < [ \k(x,y) - k(£,y)\ \f(y)\ dy < sup \k(x,y) - k(£,y)\ \\f\W . 
J I y ei 

Due to the uniform continuity of k in / x I, we have 

lim sup \k(x, y) ~ fc(£, y)\ = , 

£^ir y£l 

from which the continuity of Kf follows. □ 

Proposition 3.2. An operator K of the form (3.1) is compact from C(I) or L 1 (7) to 
either of the two spaces. 

Proof. Follow the proof of [5, Thm. 2.10] (or [15, XVII. 4]), where this is shown for 
L 2 (7) C L 1 (7), and use Holder's inequality whenever the Cauchy-Schwarz inequality is 
used. Alternatively, see [20, Exm. 3 in Sec. IV. 10]. □ 

Thus, if in our case the functions r and u are continuous, also the kernel k a is, for every 
a > 0. It then follows from Proposition 3.1 that the equilibrium density p is continuous 
as well. Therefore we can restrict our attention to C(J) in our quest for a solution of the 
eigenvalue equation (2.2). This makes the Nystrom method applicable as a discretization 
procedure, which will be presented now. 

3.1. The Nystrom method 

The Nystrom method is based on quadratures, which are used for numerical integra- 
tion, cf, e.g., Kress [14, Ch. 12]. We will use this (slightly restricted) 

Definition 3.3. A quadrature rule Q n is a mapping of the form 

Q n : C(J) - R , / Q n f = J2 a n ,kf(t n ,k) , 

fc=l 

with n e N, N n 6 N, quadrature points t n ^ E I, and quadrature weights a n .k > 0, for 
k G N n := {1, . . . , N n }. A sequence of quadrature rules, or simply a quadrature, (Q n ) is 
said to be convergent if 

Qnf - Qf for all / G C(J), (3.2) 

where Q: C(I) — > R is the linear functional that assigns to each / its integral, i.e., 
Qf = J I f(x)dx. 

Another notion that is important for the Nystrom method is the collectively compact 
convergence of operators. The standard reference for this matter is [1]. 

Definition 3.4. A sequence (K n ) of (compact) operators in a Banach space X is 
collectively compact if the set {K n B : n G N} is relatively compact (i.e., its closure is 
compact) for every bounded set B C X. If furthermore the sequence converges pointwise 
to an operator K one speaks of collectively compact convergence, in symbols K n ^> K. 
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As a direct consequence of this definition, K is compact (as well as all K n ). The central 
result for the Nystrom method is 

Theorem 3.5. Let K be a compact kernel operator of the form (3.1) whose eigenvalue 
equation 

(K - v)g = (3.3) 

is to be approximated. To this end, let (Q„) ngN be a convergent quadrature with the 
notation as in Dehnition 3.3. A complete discretization is given by the N n x N n matrices 
K n with entries 

Kn,k£ — Ci n J k(t n ,k,tn,e) , 

a partial discretization by means of the operators K n on C(I) with 

N n 

(K n f)(x) = J2 On,fc Kx, t„,fc)/(t„,fc) = Qn{k(x, .)/) . (3.4) 

k=l 

Consider the corresponding eigenvalue equations 

(K n - v n )g n = and (K n - v n )g n = , (3.5) 

where g n is an N n -dimensional vector with components g n ,k, and g n G C(I). Then, under 
the above conditions the following statements are true: 

(a) Both eigenvalue equations in (3.5) are equivalent and connected via 

9n(x) = ^2a n , k k(x,t ntk )g n , k . (3.6) 
fc=i 

(b) For every v ^ from (3.3) there is a sequence (v n ) of eigenvalues of (3.5) such 
that v n — > v as n — > oo. Conversely, every non-zero limit point of any sequence 
{v n ) of eigenvalues of (3.5) is an eigenvalue of (3.3). 

(c) Every bounded sequence (g n ) of eigenfunctions of (3.5) associated with eigenvalues 
v n —* v 7^ contains a convergent subsequence; the limit of any convergent subse- 
quence (g ni )i is an eigenfunction of (3.3) associated with the eigenvalue v (unless 
the limit is zero). 

Proof, (a) is the statement of [14, Thm. 12.7] or [5, Lemma 3.15]. (b) and (c) rely on 
K n K, which is shown, e.g., in [1, Props. 2.1, 2.2], [14, Thm. 12.8], or [5, Thm. 3.22]. 
The statements then follow from [1, Thms. 4.11, 4.17]. □ 

We will restrict ourselves to quadratures that allow for disjoint partitions of I with 
intervals I n>k , i.e., I n ,k^I n ,e $ an d Uj£o ^n,k = I, such that t n _ k e I n>k and |7„ ife | = a n ^ k 
(with k G M n )- For such quadratures it is easy to see that* 

\\Qn\\=Y, a ^ = \I\ (3-7) 
k=l 

* If not noted otherwise, the following convention for operator norms is used. If an operator maps a 
space X into itself, we denote its norm by the same symbol as the norm of X, e.g., ||.|| x , or for L 1 ; 
in all other cases the unornamented symbol |.| is used. 
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and that the partitions are unique (up to the boundary points of the intervals). Further- 
more we have 

Lemma 3.6. Let (Q n ) be a convergent quadrature that allows for partitions of I as 
described above. Then lim^oo maxj.^ \I n ,k\ = 0. 

Proof. Assume the contrary. Then there are an e > and sequences (n i ) i and (k i ) i 
with linij^oo rii = oo such that |/ ni ,fc;| > e. Due to the compactness of /, these can 
be chosen in a way that lim^oo t nit ki =: t exists. Now consider the continuous func- 
tion f(x) = max{l — 2\x — t\/e, 0}. For this we have Qf < e/2, but lim^oo Q ni f > 
elim^oo f(t n ^ki) = £, which contradicts the convergence of the quadrature (3.2). □ 

3.2. Application to the COA model 

In our case of the COA model with a compact interval / and continuous functions r 
and u, the complete discretization is given by the following N n x N n matrices: 



A n = T n -U n , K an = U n (T n -fa) 1 for a > - min w(t n fc ). 

feeA/„ 

The eigenvalue equations to be solved are 



Here, — A n + c is positive with a suitable constant c. We further have to assume that the 
A n are irreducible (which might not be the case for special choices of the t ni k, e.g., if 
ui{t n ,k) = for some k). Then, due to the Perron-Frobenius theorem, there exist (up to 
normalization) unique positive p n belonging to the eigenvalues — A„ = — p{— A n +c) +c, 
where p{M) denotes the spectral radius of a matrix M. With q n = (T n + \ n )p n also 
the eigenvalue equations 



are solved (and vice versa), cf. Lemma 2.1. 

Both K\ n:n and q n can be embedded into C(/) as described by (3.4) and (3.6). Then, 
with Theorem 3.5, one might conclude the convergence \\q n — qW^ — > 0. In the end, 
however, we are interested in the population vectors p n and their convergence to the 
density p. It might be easiest to interpret the vectors p n as point measures on /. But 
then the best one can hope for is weak convergence since the set of point measures 
is closed under the total variation norm. It will turn out that we can indeed achieve 
norm convergence if we embed the p n into L 1 (J) the following way. We choose a disjoint 
partition of I as above and let 



T n ,u = 5kt w{t n ,k) > , 
U n ,ki — a„,f u(t ni k,t n ,e) > , 



(3.8) 
(3.9) 



(A n + K)Pn = 



with p n > 0. 



(K Xn , n - l)q n = 



(3.10) 




fe=l 
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where lj denotes the characteristic function of a set J. (Note that p n k denotes the fc-th 
component of p n £ M. Nn , whereas p n is an L 1 function.) Thus the p n can be interpreted 
as probability densities on /, if we normalize them such that \\p n \\i — 1. This is most 
easily expressed using the induced norm ||/||/ n ) := J2k=i a n,k\fk\ on R Nn . Convergence 
in total variation then corresponds to \\p n — p\\ l — > [18, Thm. 6.13].* 



3.3. Convergence of eigenvalues and eigenvectors 

We now come to prove the main approximation result: 

Theorem 3.7. With the notation and assumptions from Sections 2 and 3.2, 

(a) lim^oo A„ = A > and 

(b) linin^oo \\p n —p\\i = 0, i.e., the probability measures corresponding to these densi- 
ties converge in total variation. 

The idea of the proof is as follows. In the following two lemmas, we first determine an 
upper and a lower bound for the A„ and conclude that there is a convergent subsequence. 
Then we show that every convergent subsequence converges to A and hence the sequence 
itself. By Theorem 3.5, this implies the convergence of a subsequence of (<Z„/||<Z„||oo) 
to a (non-negative) limit function. Since, due to Theorem 2.2, the latter is unique, we 
conclude that it is ^/H^Hoo- With this, part (b) can be shown. 

Lemma 3.8. There is a constant M > such that |A„| < M for all n E N. 
Proof. Using (3.8) and (3.9), one checks 



l\iPnll(n) \\-A- n P. 



. U"nrnU(n) » nfn II in) _ \ ^ 

K\ = -n — n = -n — < su p 2^ 1 

HPnll(n) WPn\\(n) ll/ll(„) =1 fc=l 



N„ 



--i 

< maxM)(t k ) +maxu(£ k ,t e ) Va k 

k k,t ' 

fe=l 

< Iklloo + hWcdxi) s up ||Q J| =: m > . 



IC(/x/) • 

Hi 

Here, ||Q m || = |/| due to (3.7). More generally, sup m ||Q m || < oo holds for any convergent 
quadrature according to the theorem of Banach-Steinhaus, compare [19, Thm. 2.5]. □ 

Lemma 3.9. liminf A„ > 0. 

71— >00 

Proof. We start by following Burger [4, p. 134] and show that the spectral radius 
p(K a ) is larger than 1 for sufficiently small a > 0, from which then A > a > follows 
by Lemma 2.3. Let J be the interval from (2.9). Then we have 

(K a lj)(x)= [ U ^ dy>\j(x)^mi u(x',y') ( (w(y) + a)' 1 dy 
Jjw{y)+a x>, y >eJ Jj 

* One may also define operator analogs of the A n , see [17, Sec. II. 2. 1.2]. 
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and thus 



||if Q m Hi 1/m > easMu(x,y) [ (w(y) + a)- 1 dy for all m G N, 

which implies for the spectral radius 

p{K a ) > essinf u(x,y) / (w(y) + a)^ 1 dy. (3-11) 
x,ye.J J j 

The RHS IS - clS cL function of a, strictly decreasing. Thus, as a consequence of B. Levi's 
monotone convergence theorem [9, Thm. III. 12.22], also 



lim p{K a ) > essinf u(x, y) I (w(y)) 1 dy 

a\0 x,ye.J J j 



> 1 



according to (2.9) (including divergence of both sides). 

Now we choose a > such that the RHS of (3.11) is larger than or equal to 1 + e, 
with a sufficiently small e > 0. Furthermore, we pick, according to the convergence of 
the quadrature, an n with ess inf^gj u(x,y)\Q n (w + a)^ 1 — Q(w + a) _1 | < e/2 for all 
n > no- This way 

(K a _ n lj)(x) = Q n (u(x,.)(w + ay 1 l J ) > lj(x) essinf u(x',y')Q n {w + a)^ 1 

x' ,y'e.J 

>lj(x) ^essi|f«(^,y , )Q(«; + a)- 1 -0 >lj(x) (l + |) . 

Hence, by Lemma 2.3, A„ > a > for all n > n , from which the claim follows. □ 

Proof of Theorem 3.7. By Lemmas 3.8 and 3.9, the sequence (A„)„ has a convergent 
subsequence (A„J i with limit A' G ]0, M]. Consider (Kyf)(x) = Q(k\>(x, .)/) as well as 
(K n f)(x) := {K Xntn f)(x) = Q n (T + X n )-\T + \')(k y (x, .)/). We first show that the 
'distorted' quadrature Q ni = Q ni (T + A„J _1 (T + A') is convergent. Note that, for i 
large enough, such that infj> io A^ > 0, and i > i , 



|(T + A„,)- 1 (T + A)-1|| 00 = sup 

ll/lloc^ 1 



" +A -/-/ 



w + \ n% 



(3.12) 



<||( w +inf A^) -1 !!^ |A - A n J ll/H -> 0. 

3>io 

Then, since (Q n ) is convergent by assumption, we have, for all / G C(I), 

II - <3)/IL < IIQ ni ((T + AnJ-^T + A') - 1X/IL + \\(Q ni - Q)/!^ - 0, 

where the first term vanishes in the limit due to sup m ||Q m || < 00 and (3.12). 

With this it follows from Theorem 3.5 that p{K n ) = 1 is also an eigenvalue of Ky 
going with a non-negative eigenfunction. The latter is even a.e. positive since, due to 
the irreducibility (2.8) of K\i, there cannot be a set with positive measure on which 
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a non-negative eigcnfunction vanishes.* But since, according to Theorem 2.2, there is, 
up to normalization, only one positive eigcnfunction, wc have A' = A. Therefore every 
convergent subsequence of (A„)„ converges to A, and thus, due to the boundedness, also 
the sequence itself. This proves part (a). 

Along the same line of reasoning, (a n q n ), with a n = l/UgJIoo, has a convergent subse- 
quence with an a.e. positive limit function, which equals aq with a = 1/ 1| oo - Therefore, 
\\ a n 1n ~ a l\\oo for n ^ oo. Now consider 

\\ a nPn - Hloo = \\ a nPn ~ ( T + ^"Vlloo 

= max sup \(w(t n k ) + A„) _1 a„g„ Jc - (w(x) + \)~ 1 aq(x)\ 
< max|(w(i„, fc ) + AJ" 1 - (w(t n k ) + A) -1 1 a n q n (t n ^ k ) + 
max (w(t k ) + \)~ 1 \a n q n (t n ^ k ) - aq(t fe )| + 
max sup \(w(t nk ) + \y 1 aq{t n k ) - (w(x) + A) _1 ag(a;)| . 

The first term is bounded from above by 

\\-\ n \\\(w+ irif A m )- 1 ( w + A)- 1 || co , 

m>no 

for n > riQ with sufficiently large no, and vanishes for n — > oo due to A„ — > A. The second 
term vanishes due to the uniform convergence of the a n q n towards aq, and the third due 
to the uniform continuity of (w + X)~ 1 q and Lemma 3.6. With this, a n p n — > ap in L°°(/) 
and thus in L 1 (7). Hence, a n — > a and p n — > p in L 1 (7), which proves part (b). □ 

4. Discretization — unbounded interval 

Now we assume the types to be taken from I = K and the functions r and u to be 
continuous. It will be one aim of this section to analyze what further conditions have 
to be imposed in order to allow for a discretization procedure similar to the one in the 
previous section. In order to do so, we start by a summary of the relevant theory. 

4.1. The Galerkin method 

In the Galerkin method, an approximation of compact operators is achieved using 
projections to finite-dimensional subspaces. This method has been reviewed, e.g., by 
Krasnoscl'skii et al. [13, Sec. 18]. The results needed in the sequel are collected in 

Theorem 4.1. Let K be a compact linear operator on the Banach space Y. Consider 
the eigenvalue equation 

{K-v)g = Q, (4.1) 

* Let q be the eigcnfunction and J = {x : q(x) > 0} with < \ J\. Assume \J\ < \I\. Then, for x g I\J, 
we have = q(x) = j j k\i{x,y)q(y) Ay, which implies, for a.e. y € J, that ky(x,y)q(y) = and thus 
u(x,y) = 0, contradicting (2.8). 
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which is to be approximated. To this end, let (Y n ) be a sequence of closed subspaces of 
Y with bounded projections P n onto them. On these subspaces, let the compact linear 
operators K n be dehned, together with the eigenvalue equations 

(K n - u n ) 9n = . (4.2) 

Assume that 

\\K n - P n K\\ Yn , \\K- P n K\\ Y as rw oo. (4.3) 
Then the following statements are true: 

(a) For every v ^ from (4.1) there is a sequence (v n ) of eigenvalues of (4.2) such 
that v n — > v as n — > oo. Conversely, every non-zero limit point of any sequence 
(u n ) of eigenvalues of (4.2) is an eigenvalue of (4.1). 

(b) Every bounded sequence (g n ) of eigenvectors of (4.2) associated with eigenvalues 
v, n — > v 7^ contains a convergent subsequence; the limit of any convergent sub- 
sequence {g ni )i is an eigenvector of (4.1) associated with the eigenvalue v (unless 
the limit is zero). 

Proof. Sec [13, Thms. 18.1, 18.2]. □ 
A sufficient condition for the validity of the second assumption in (4.3) is given by 

Proposition 4.2. Let X be a normed space, Y a Banach space, and K: X — > Y a 
compact linear operator. For bounded linear operators P n : Y — > Y (n 6 N) with P n — ► 1 
pointwise for n — > oo, the operators P n K approximate K, i.e., \\P n K — K\\ — > 0. 

Proof. Follow the proof of [21, Thm. II.3.5], where the additional assumptions on X 
and (P n ) are not used. □ 

4.2. Application to kernel operators 

In our case of the COA model we have X = Y = L 1 (K) and K is of the form 

(Kf)(x) = f k(x, y)f(y) dy for all x e R (4.4) 

with a measurable kernel i:MxM^I. Therefore, for the Galerkin method to work, it 
is necessary that, for L 1 (K), operators P n as in Proposition 4.2 exist. We will explicitly 
construct such operators using a sequence ({I n ,k ■ 1 < k < N n }) n of families of disjoint 
intervals that get finer and finer and also ultimately cover every bounded interval.* 

Proposition 4.3. Let Y be the Banach space L 1 (R) and finite- dimensional subspaces 
Y n of Y chosen to consist of all step functions with prescribed (bounded) intervals I n ,k 
(k G J\f n := {1, . . . , N n }) with the following properties: 

* Both properties are formally captured by (II) in Proposition 4.3. 
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(11) For every bounded interval I C R and every e > there is an no such that, for all 
n > n , a set L C N n exists for which I n L := {J eeL I n ,e satisfies \I\I n ,L\ — and 
\Iu,l\I\ < £■ (We then say that I is e-optimally covered.) 

(12) |7„ )fc n I n>e \ = for all n eN and 1 < k < £ < N n . 

Then, with the characteristic functions <~p n ^ = 1 j , tic projections P n onto the sub- 
spaces Y n spanned by {ip n ,k ■ k G 7V„} are given by 

where Jj ^ f(x) dx are conditional expectations (compare [20, Thm. IV. 2. 4]*). The pro- 
jections satisfy \\P n \\i = 1 an d P n — > 1 pointwise. 

Proof. Obviously, the subspaces V„ are closed, finite-dimensional, and the P n are, 
due to (12), projections onto them. Since 

\\p n f\\ = £ / 1/(^)1 da; < / da; = ll/lli fOT ever y / e Ll ( R ) 

fc=l ■'•fn,* ■'H 

and ||P„v„,fe|| = ||<£n,fc|| for every fc G 7V„, we have ||P„|| = 1. 

We now show that P n — > 1 pointwise. To this end, let / G L X (K) and e > be given. 
Remember that the set of all step functions is, by definition, dense in L 1 (M), compare 
[15, Sec. VI. 3]. Therefore, we can find a step function tp = J2h=i ^ k ^-J k ( Wl th bounded 
intervals J&) that satisfies ||/ — V'lli < Due to (II) we can now choose an n such that 
| Ul!Li Jk\ Ufe^i Ai,fe| — for all n > no- Then, the only contributions to \\P n ip — are 
due to mismatches at the boundaries of the Jk- Therefore, let and (k G N n ) be 
open intervals of measure n = e/(12mmaxfc \ipk\) that contain the right and left boundary 
points of Jfe, respectively. Choosing n\ > no according to (II) large enough such that 
every is 77-optimally covered for n > n\, we have \\P n ip — tp\\i < 2^a~i ^Vpk < e/3 
for n>n\. Putting everything together yields, for n > ni, 

||P„/ - /Ik < ||P n (/ - VOIIi + - V'lll + U ~ /111 < 

which proves ||P„/ — /|| — > for n — > oo and thus the approximation property. □ 

With respect to a kernel operator if of the form (4.4) and some / = J2k=i Vn,kSk in 
Y n , the above procedure amounts to the discretization 

[P n Kf) = Y]^n,kT^— | f / / k(x,y)Y]f e ip n j(y)dydx) 

k=1 \-ln,k\ \Jl nik JR t=1 ) 



k=i 



1yf—\[ I I k(x,y)dydx J ,/ £ =: ^ ¥Vfc Z M n.ktfi 

i |J ™> fc V- 7 '".*.- / fe=i £=1 



* See also [17] for a discussion of the connection to the approximation property of Banach spaces. 
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with an N n x N n matrix M n — {M n ^t)- The corresponding eigenvalue equation is 

M n9 n = v n9n ! or cquivalently P„i^„ = v n g n , 

where j„ e F n is granted due to the projection property. An example of intervals I n ,k sat- 
isfying (II) and (12) is J n , fe = [-n+2" n (fc-l), -n+2~ n k] with k G A/"„ = {1, . . . , 2 n+1 n}. 
With respect to compactness of K, following Jorgens [11, Sees. 11, 12], we extract 

Proposition 4.4. A kernel operator K on L X (R) of the form (4.4) is compact if it 
satisfies the following conditions: 

(CI) The function x — ► fc(x, .) from R to L X (R) is continuous and bounded. 

(C2) For every e > there exists a finite open covering (V\, . . . , V^) of R and points 
Xj G V, such that ||fc(x, .) — k(xj, < e for all x G Vj and ah" j. 

(C3) The function y — ► k(.,y) from R to L X (R) is continuous and bounded. 

(C4) For every e > there exists a finite open covering {W\, . . . , W^) of R and points 
yj G such that ||fc(.,y) — A;(.,yj)|| 1 < e for ah" j/ G and ah j. 

Proof. First, as in [11, Sec. 12.4], we consider the dual system (C(R),Ci(R)) with 
the bilinear form (/, g) = J f(x)g(x) dx. Here, C(R) is equipped with the supremum 
norm H-H^ and Ci(M) := C(R)nL 1 (R) with the norm |||.||| := max{||.|| 00 , ||-||i}. With 
this, we dehne the transposed K T oi K via (K T g)(y) — J R g(x)k(x, y) dx, for all y G R. 
Then, by (C1)-(C4) and [11, Thms. 12.2, 12.3], the compactness of K and K T on C(R) 
follows. 

As both K T and K are bounded as operators on C(R), they are, at the same time, Hille 
Tamarkin operators in 7Yoooo(R) since the respective norm, l-l^oo in (2.7), is just given by 
sup^gjj J R k(x, y) dy and sup yCR / R k(x, y) dx, respectively [11, Thms. 12.2, 12.3]. Then, 
according to [11, Thm. 11.5], K and K T can also be regarded as bounded operators on 
L X (R); thus, both map Ci(R) into itself. Due to [11, Thm. 12.6] there is, for every e > 0, 
an operator of hnitc rank, K e , with |||.K" e — K\j < e, where |||A||| := max{||A|| 0O , ||A T || 00 } 
is a norm for the Banach algebra of all operators on C(R) that map Ci(R) into itself 
and have a transposed of the same kind. We have |||A/||| < |||A||||||/||| for / G Ci(R), see 
[11, Sec. 12.4]. Thus, ||| A||| can serve as an upper bound for the operator norm of A on 
Ci(R). Therefore, K is compact as an operator on Ci(R) and can be approximated by 
K e . Furthermore, according to [11, Thm. 11.5], \\K E - < \\(K e - K) T \\ < e holds. 
Hence, K is compact as an operator on L X (R) as well. □ 

4.3. Application to the COA model 

Checking the compactness of K a by conditions (C1)-(C4) of Proposition 4.4, we would 
be able to apply Theorem 4.1 and approximate K a by operators of hnitc rank. However, 
the original system is described by the (non-compact) operator A = T — U , not by some 
K a . It will be shown that it is indeed possible to discrctize the operators T and U directly 
by applying the projections P n from Proposition 4.3, if further restrictions apply. Then, 
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even more generally, the approximation can be done by choosing arbitrary points in the 
intervals I nk at which the functions w and u are sampled. Both procedures will now be 
described in detail. 

In the first setting, K\ is approximated by K n :— P n U(P n T + A„) _1 . Explicitly, for 
/ G Y n with / = Ylk=i fk<Pn,k, it reads 

N n 1 N n 

P n Tf= Vn.fcJT 1 / w(x)dxf k = Y fnMw(t^j,)f k , 

fe=l |J "' fc| Jl n,k fe=l 

k,l=l \ ±n M JI n ,k JI n ,e M=1 
with appropriate points t™ k , t™ kl G I nk and t^ y ke G /„ £ that satisfy 

Tj—; I <x, tl%) Ax = u(t^ M , tl%) . (4.5) 

These exist due to the continuity of w and u. But more generally, we may pick the points 
arbitrarily from the respective intervals. 

In either case, we define the N n x N n matrices T„, U n , and A n := T n — U n via 

T n ,kk := w(tZ k ) , U nM := \I n A u{t u n %X%) ■ (4-6) 
The corresponding operators in Y n are given by 

Tnf = '^2 t Pn,kT n ,kk.fk , U n f = ^ ^fn,kU n ^klft , A n = T n — U n , 

k=l k,i=l 

again with / = Y^k=i fk<Pn,k- For notational convenience, we also define the matrices 
P a ,n by 

PnK a f = Y] ipn.k-rr-^ [ [ U ^ X,y ^ dydxfi-. V tpn.kPa.n.ktfa ■ (4.7) 

^ \I n ,k\ Ji n ,Ji n , e w(y)+a ^ 
The eigenvalue equation to be solved is 

(A n + \ n ) Pn = 0, 

which is equivalent to 

(A n + X n )p n = 0, (4.8) 

where p n = Pn,fc</>n,fc e Y n- with K a,n = U n (T n + a)- 1 , a > -mm keAfn w(t n , k ), 

and q n = (T n + \ n )p n also the eigenvalue equation 

(K Xn , n - l)q n = 
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is solved (and vice versa), cf. Lemma 2.1. (The inequality A„ > — mmkeJ\r n w(t n: k) follows 
from Theorem 2.2.) 

For these procedures to be valid approximations, the first condition in (4.3), that is, 
\\K n — PnK\\y ~* 0> nas to be true for K = K\ and K n = C/„(T„ + A„) _1 . This, however, 
is not given automatically. Problems arise from the fact that in K n the averaging defined 
by P n (or, more generally, the sampling) is applied to the enumerator and denominator 
of k\ n separately, whereas in P n K the quotient k\ is averaged as such. It turns out that 
some additional requirements of uniform continuity are sufficient for the convergence. 
This is made precise in the following two propositions. 

Proposition 4.5. Suppose that the following conditions are true: 

(51) u(x, .) is uniformly continuous for all 

(52) k a is uniformly continuous on I x R for all a > and all bounded Id. 

(53) There is a function w min : R — * R>cb satisfying 

sup yAr dx < oo for all a > 0, 

rs/gR w min (y) + a 

and an n G N such that w(y) > w m i n (y') for all n > no, £ G Af n , and y, y' G I n ,t- 

Then, for K = K a and K n = P n U(P n T + a) -1 with any a > and the projections 
P n from Proposition 4.3, the first condition in (4.3) is fulfilled, i.e., \\K n — P n K\\ Y —> 
as n — > oo. The same is true for K n = K a<n = U n (T n + a) -1 with the more general 
discretization from above if in addition to (S1)-(S3) the following condition is satisfied: 

(54) There is a function u max :MxM^ R>o? satisfying 

sup — dx < oo for all a > 0, 

y eR w min (y) + a 

and an m > n such that u(x, y) < u max (x' , y) for all n > m, k G N n , t/GE, and 

X, X G In,k- 

Let us split the rather technical proof into a couple of digestible lemmas. 

Lemma 4.6. If conditions (SI) and (S2) are true, then for every e > and every 
compact interval I C M there is an n 2 such that for all n > n 2 and all k, £ G N n with 
I n ,k H / ^ we have 



1 



+ a 



e 



Proof. Let e and I be given as above and I = U n eN Ufe-/„ k ni^Hi which i s a 
bounded interval due to (II) from Proposition 4.3. By assumptions (SI) and (S2), u and 
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k a are uniformly continuous on Iq x R. Further, (w + a) -1 is bounded by oT 1 . Thus, 
there is an n 2 such that, for every n > n 2 and k, I G M n with I n k 0/^0, 



\I n ,k\ \I n ,e\ Ji n , k Ji n>e w(y)+a w{t™ t ) + a 



< 



Here, the points t£ x fe< G 7„ ; fe and G are chosen such that the first equality holds, 
which is possible due to the continuity of k a . From this the claim follows easily with (4.6) 
and (4.7). □ 

Lemma 4.7. For every e > there is a compact interval I\ such that, for all intervals 
I D h and all n G N, 



i 1 ct.n.Kt 

k max — — < e . 



*6AT„ |/ n ,<| 



/„. fc n/=0 



Proof. Due to (U4) there is a compact interval I\ such that, for all I D I\, 

p 

1 a.n 



V |J», fc | max P ^ n ' ke < V" |j fe |_!_ / 
^ 'eeXn \I n , e \ k ' ™' fe l 

j n , fe n/=0 / n , fe n7=0 



< 



/ 



max 



u{x,y) 
max — — ax 

yevi w(y) + a 

u(x,y) 



RVl ym w{y) + a 



dx < e . 



which proves the claim. 



□ 



Lemma 4.8. If condition (S3) is true, and if 

(i) U nM = J In k u(x, t u n v ke ) dx for all k,£eAf n or 

(ii) condition (S4) is fulfilled, 

then for every e > there is a compact interval I2 such that, for all intervals I D I 2 and 
all n G N, 



El t I U n kl 
\l n fe max — — — 
; 1 n ' fc| /&V„ |J„,/|(T n ,« + a 

J„,*n/=0 



< e. 
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Proof. In case (i) we have, using (4.5), 

^ , T , x ^ Cm) ds 
> L i. max — — — - — ■ — = > max : — ; 

< y f max u fo g > dx < / max ^ dx < e 



, , k yGR W min (y) + a ' 7 R \ 72 yGR W min (y) + OL 

i n , k m=® 

for some compact interval I2, due to (S3), and all intervals / D li- In case (ii) we can 
find, due to (S4), a compact interval 7 2 such that, for all intervals I D I2, 

\L„ h max < > i„ i. max- 



teN n w(t™j)+a ^ ' yGR w m in(y) + a 

/n,fen/=0 /„, fc n/=0 

N //„ ,„ u max(a;, y) da; ^ ^ u max (x,y) 



max ; — < / max — ax < e . 

k yGR w min (y) + a J R \ l2 yGR w min (y) + a 

i n , k m=® 

This completes the proof. □ 

Proof of Proposition 4.5. Let e > be given. Choose a compact interval / such 
that I D Ii U I2 with I\ and I2 from Lemmas 4.7 and 4.8. Let I3 ~ [J n k .j kr] i^$ In.k- 
Further, let no be as in (S3), n\ as in (S4) (or no = ri\ if not applicable), n 2 as in 
Lemma 4.6, and n > max{n , ni, n 2 }. Then 



(Pnifa - -£fa,n||y n = Sup y \I n 



f£Y r 



fc=l 



Pa,n,fe« - 7^ L " ft 

Tn.ee + a ' 



1=1 



< 



( ? + ? ) |J - fc| ^]i 



^ l^nfelv— t+ l-fra fel I max — — — ^ + max , "' fc£ ) < 3e 

i n , fe n7#0 /„, fc n/=0 

according to Lemmas 4.6-4.8. From this the claim follows. □ 

Proposition 4.9. Let a n > — minfee^ w{tn,k) w ^ n a n — >a>0asn— >oo and the 
hypotheses of Proposition 4.5 be satisfied. Then \\K an ^ n — P n K a \\ Yn — > 0. 

Proof. Consider 



\\P n K a -U n (T n + a n y 1 \\ Yn 

< \\P n K a - UniTn + a)- 1 ^ + ||C/„[(T„ + a„)- 1 - (T„ + a)" 1 ] l|y„ 
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The first term tends to zero as n — > oo according to Proposition 4.5. For the second, 
choose n such that inf„>„ a n > 0. Then, for n > n , 

WU^ + a^ 1 -{T n + a)- 1 ]\\ Yn 

= \a-a n \\\U n (T n + a n )- 1 (T n + a)- 1 \\ Yn < \a - a n \ \\U\\ Y { inf aj- 1 ^ 1 . 

n>no 

This vanishes as n — > oo since all constants that occur are finite, from which the claim 
follows. □ 

4.4. Convergence of eigenvalues and eigenvectors 

Let us now show 

Theorem 4.10. With the notation and assumptions from Section 2 and A, p, X n , p n 
as in (2.2) and (4.8), we have 

(a) lim„^oo A„ = A > and 

(b) lirn„_ >00 \\p n —p\\i = 0, i.e., the probability measures corresponding to these densi- 
ties converge in total variation. 

The plan is the same as described in Section 3.3. The proofs, however, are quite different 
due to the more general setup. 

Lemma 4.11. There is a constant M > such that limsup^^ A„ < M. 

Proof. Choose an a > such that ||if Q || y < 1 — e for some < e < 1, which is 
possible since ||/r Q ,||y — > for a — > oo. Then, for all n > no with some no, due to 
Propositions 4.2 and 4.5, |||P n Jf a ||y-||Jf a || y | < e/3 and |||#a,„||y B - \\P n K a \\ Yn \ < e/3. 
For these n, we have 

p(K a , n ) < \\K a , n \\ Yn < \\P n K a \\ Yn +e/3 < \\P n K a \\ Y + e/3 

< \\K a \\ Y + 2s/3< l-e/3<l 

and thus A„ < a by Lemma 2.3. Then, with M = a, the claim follows. □ 

Lemma 4.12. liminfn^ooAn > 0. 

Proof. In a modification of the proof of Lemma 3.9, we choose a > such that 
p{K a ) > 1 + £ with a sufficiently small e > 0. We know from the theorem of Jentzsch 
[20, Thm. V.6.6] that p{K a ) is a simple eigenvalue of K a and the only one with a positive 
eigenfunction. The same is true for p(K a>n ) with respect to K a<n (as an operator in Y n ). 
Theorem 4.1 together with Proposition 4.5 implies that there is a sequence of eigenvalues 
y n of K aiTl with limit p(K a ). Therefore, liminfn^oo p(K a-n ) > p(K a ) > 1 + e and thus 
A„ > a > for sufficiently large n. From this the claim follows. □ 
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Proof of Theorem 4.10. From Lemma 4.11 and 4.12 we conclude that there is 
a convergent subsequence (A„J i with limit A' G ]0, M]. Then, due to Proposition 4.9, 
K\ nv m converges to P n K\> in norm. Hence, lim^x, p(K\ n . ) = p(K\ ni>ni ) = 1 is 
an eigenvalue of K\* by Theorem 4.1. Furthermore, a subsequence of (a n .q n .), where 
a n = l/||9„||i, converges to an eigenfunction q of Ky, and q > (but g^O). As there 
is only one non- negative eigenfunction by Theorem 2.2, we conclude A' = A and q = aq 
with a = 1/Hgllp Since this is true for every convergent subsequence of (A„), the claim 
of part (a) and the convergence a n q n — > aq follow. 

Now, let no be sufficiently large such that a := inf„>„ A„ > 0. Then, for n > no, 

\KPn -(T + A)- 1 ag|| 1 = ||(P„T + X n )~ 1 a n q n - (T + A)" 1 ^^ 

<||[(P n T + A„)- 1 -(P n T + A)- 1 ]a n gJ| 1 + 

||(P n T + A)- 1 (a„ <Z „- a(Z )|| 1 + 

H^T + A^-CT + A)- 1 ]^! 
< ^xlA - AJ + {\\a n q n - aq^ + £||(1 - P n) Ta 9lli - °- 

With this, a„p„ — > ap in L 1 (7), hence a„ — > a and p„ — ► p, which proves part (b). □ 
5. Comparison of both methods 

Both approaches, the application of the Nystrom method in the case of a compact interval 
and of the Galcrkin method in the case of an unbounded interval, effectively lead to 
the same approximation procedure in our case of the CO A model. First, one chooses 
appropriate intervals l n .k and points t nt k G I n ,k (also for an unbounded interval the use 
of identical points t™ k = t™ x ke = t^ ek = t n ^ k seems reasonable in many cases). Then, 
the operators T and U from (2.3) and (2.4), respectively, are approximated by matrices 
T n and U n , cf. (3.8), (3.9), and (4.6). For these, the (finite-dimensional) eigenvalue 
problem (T n — U n + \ n )p n = is solved. Here, the eigenvectors p n are considered as 
probability densities on /. Then, under the conditions described above, the eigenvalues 
A„ converge to A and the measures corresponding to the p n converge in total variation 
to the equilibrium distribution described by the solution p of the original problem (1.1). 

The differences between the two approaches lie on the intermediate technical level of the 
compact operators K a and K a ^ n and the solutions q and q n of the equivalent eigenvalue 
problems (2.6), (3.5), and (4.2). Here, in the first case we have collectively compact 
convergence K\ n . n ^> K\ going together with \\q n — qW^ — > 0, whereas in the second 
case \\P„Kx - K X \\ Y -> in Y = L X (R) and \\K Xn ,„ - P n K x \\ Yn -> in the subspaces 
Y n going together with \\q n — q\\ ± — > 0. On this level, neither does ||^A„,n — ^a||oo — ^ 
hold in the first case, compare [14, Thm. 12.8], nor any kind of collectively compact 
convergence in the second. 

Both methods may, strictly speaking, only be applied to continuous mutation ker- 
nels u. This excludes, for example, T-distributions (reflected at the source type), where 
u(x,y) cx \x — y| e_1 cxp(— d \x — y\), which have poles for x — y if 9 G ]0, 1[ and 
d > 0. These distributions incorporate biologically desirable properties, such as strong 
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leptokurticity, and have been used, e.g., in [10]. However, kernels as the above may be 
approximated arbitrarily well by continuous ones in the sense that the norm of the differ- 
ence operator — and thus the difference of the largest eigenvalues — gets arbitrarily small. 
Then, the procedures described here may be applied to these continuous kernels. 

6. Outlook 

This article shows that most reasonable COA models can be approximated arbitrarily 
well by models with discrete types. Therefore, one can expect both model classes to 
behave quite similarly. For certain mutation-selection models with discrete types, a sim- 
ple maximum principle for the equilibrium mean fitness A was recently found [8] (see 
also [7, 6, 2]). It takes the form 

A ~ sup(r(x) — g(x)) 
xei 

and holds as an exact identity in a limit of infinitely many types that densely fill a compact 
interval /. In the simplest linear ordering of types is assumed and mutation is 

taken to only connect every type x with its two neighbors at rates u ± (x). Then, the 
function g is given as g(x) = u + (x) + u~(x) — 2yju + (x) u~(x). In a subsequent analysis 
[6] , models with three types of mutation — and hence six neighbors of every type — were 
considered. For these, g is given as the sum of three terms of the above pattern (and x 
has three components), one for each type of mutation. 

In the light of the findings presented here, one may conjecture that also for certain 
COA models the above characterization is valid with an appropriate function g. First 
steps in [17], both analytical and numerical, corroborate this conjecture with 

9{x) = J (u(x, y) - yju(x,y)u(y,x)) dy , 

which generalizes the additive structure of g found in [6] with respect to a continuum of 
possible mutations. The important prerequisite seems to be the possibility to approximate 
every local subsystem, corresponding to a small interval J C I, by a COA model whose 
mutation kernel is of the form u(x, y) = cxp(7 (x — y)) h(\x — y\). Then, in a limit v — > oo, 
where h is replaced by h v (\x — y\) = v h(y \x — y\), the above expression seems to become 
exact. A rigorous proof for this statement seems feasible in the near future. 
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